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Spin-glass phases and phase transitions for q-state clock models and their q — ^ oo limit the XY 
model, in spatial dimension d — 3, are studied by a detailed renormalization-group study that is 
exact for the d — 3 hierarchical lattice and approximate for the cubic lattice. In addition to the now 
well-established chaotic rescaling behavior of the spin-glass phase, each of the two types of spin-glass 
phase boundaries displays, under renormalization-group trajectories, their own distinctive chaotic 
behavior. These chaotic renormalization-group trajectories subdivide into two categories, namely as 
strong-coupling chaos (in the spin-glass phase and, distinctly, on the spinglass-ferromagnetic phase 
boundary) and as intermediate-coupling chaos (on the spinglass-paramagnetic phase boundary). We 
thus characterize each different phase and phase boundary exhibiting chaos by its distinct Lyapunov 
exponent, which we calculate. We show that, under renormalization-group, chaotic trajectories 
and fixed distributions are mechanistically and quantitatively equivalent. The phase diagrams of 
arbitrary even q-state clock spin-glass models in d = 3 are calculated. These models, for all non- 
infinite q, have a finite-temperature spin-glass phase. Furthermore, the spin-glass phases exhibit a 
universal ordering behavior, independent of q. The spin-glass phases and the spinglass-paramagnetic 
phase boundaries exhibit universal fixed distributions, chaotic trajectories and Lyapunov exponents. 
In the XY model limit, our calculations indicate a zero-temperature spin-glass phase. 
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I. INTRODUCTION 

Spin-glass phases, with randomly frozen local order [l[ 
and chaotic behavior under scale change Q , reflecting the 
effects of frozen interaction disorder, competition, and 
frustration, remain a uniquely fascinating and broadly 
relevant subject of statistical mechanics and condensed 
matter physics. However, the large and richly complex 
amount of theoretical knowledge produced on spin glasses 
has been overwhelmingly derived from Ising, i.e., Sj = 
±1, spin models. Q 

By contrast, we present here a detailed renormaliza- 
tion-group study of spin-glass phases and phase tran- 
sitions, for g-state clock models and their q — > oo 
limit the XY model, in spatial dimension d = 3. 
We note that, in addition to the now well-established 
chaotic behavior of the spin-glass phase @, I34l6|. each 
of the two types of spin-glass phase boundaries dis- 
plays, under renormalization-group trajectories, their 
own distinctive chaotic behavior. We see that these 
chaotic renormalization-group trajectories subdivide into 
two categories, namely as strong-coupling chaos (in 
the spin-glass phase and, distinctly, on the spinglass- 
ferromagnetic phase boundary) and as intermediate- 
coupling chaos (on the spinglass-paramagnetic phase 
boundary) . We thus quantitatively characterize each dif- 
ferent phase and phase boundary exhibiting chaos by its 
distinct Lyapunov exponent, which we calculate. We 
show that, under renormalization-group, chaotic trajec- 
tories and fixed distributions are mechanistically and 
quantitatively equivalent. 

We calculate and display the phase diagrams of arbi- 
trary even g-state clock spin-glass models in d = 3. These 
models, for any non-infinite q, have a finite-temperature 



spin-glass phase. Furthermore, we find that the spin- 
glass phases exhibit a universal ordering behavior, inde- 
pendent of q. The spin-glass phases and the spinglass- 
paramagnetic phase boundary exhibit universal fixed dis- 
tributions, chaotic trajectories and Lyapunov exponents. 
In the d = 3 XY model limit, our calculations indicate a 
zero-temperature spin-glass phase. 



II. THE q-STATE CLOCK SPIN-GLASS MODEL 
AND THE RENORMALIZATION-GROUP 
METHOD 

The g-state clock models are composed of unit spins 
that are confined to a plane and that can only point along 
q angularly equidistant directions. Accordingly, the q- 
state clock spin-glass model is defined by the Hamiltonian 
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where /3 = 1/ksT, at site i the spin angle 0i takes on 
the values (2Ti/q)ai with <n — 0,1,2, ...,q — 1, and (ij) 
denotes that the sum runs over all nearest-neighbor pairs 
of sites. The bond strengths Jjj are + J > (ferromag- 
netic) with probability 1 — p and — J (antiferromagnetic) 
with probability p. This model becomes the Ising model 
for q=2 and the XY model for q — > oo. 

The g-state clock spin-glass model, in d = 3 di- 
mensions, is readily solved by a renormalization-group 
method that is approximate on the cubic lattice [171 . 1 1 811 
and simultaneously exact on the hierarchical lattice [L9l — 
[23j . Under rescaling, for q > 4, the form of the interac- 
tion as given in the rightmost side of Eq.(l) is not con- 
served and one must therefore express the Hamiltonian 
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FIG. 1: (a) Migdal-Kadanoff approximate renormalization- 
group transformation for the d = 3 cubic lattice with the 
length-rescaling factor of 6 = 3. Bond-moving is followed by 
decimation, (b) Exact renormalization-group transformation 
for the equivalent d = 3 hierarchical lattice with the length- 
rescaling factor of 6 = 3. (c) Pairwise applications of the 
quenched probability convolution of Eq.(5), leading to the 
exact transformation in (b). 



more generally, as 
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where, with no loss of generality, the maximum value of 
V(9i — 6j) is set to zero. 

The renormalization-group transformation, for spatial 
dimensions d = 3 and length rescaling factor b = 3 (nec- 
essary for treating the ferromagnetic and antiferromag- 
netic correlations on equal footing), is achieved by a se- 
quence of bond moving 
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FIG. 2: (Color online) Calculated phase diagrams of the 
g = 2,4, 6, 12 clock spin-glass models in d — 3 dimensions. 
The phase diagram for the XY limit, namely q — > oo, is also 
shown, by the dashed curve, calculated here with q = 360 
clock states. As q is increased, it is found that the spin-glass 
phase retreats to lower temperatures while further protruding 
into the ferromagnetic phase. In the XY limit, the spin-glass 
phase disappears at zero temperature (see Fig. 4 below), as 
the phase boundary between the remaining ferromagnetic and 
paramagnetic phases numerically stabilizes, on the scale of the 
figure, for q > 6. 



\r in a \ \ n m la a \ io\ where the constants Gu ; are fixed by the requirement that 

Vbm\Pl ~ 02) + Cri2 = y J V n {t>i - fl 2 J (3) . « ? 

~ J the maximum value ot V (Ui — Uj) is zero. 

and decimation 

The starting bimodal quenched probability distribu- 

e v dec (e 1 -e 4 )+G 1 4 = e Vi{0i-02)+v 2 (e 2 -e 3 )+v 3 (e 3 -e i ) ^ tion of thc mteractionSj characterized by p and described 

6 2 ,6 3 above, is also not conserved under rescaling. The renor- 

(4) malized quenched probability distribution of the interac- 
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FIG. 3: (Color online) Phase diagrams of the q — 
12, 18, 36 clock spin-glass models in d = 3 dimensions. 
Two new phenomena are simultaneously detected here: 
(1) double reentrance: paramagnetic-ferromagnetic-spinglass- 
ferromagnetic as temperature is lowered, (2) lateral reen- 
trance: ferromagnetic-spinglass- ferromagnetic-paramagnetic 
as p is increased. The XY limit is given, calculated here 
with q — 360 clock states, by the dashed curve. The panels 
on the right show the region of the multicritical point, where 
the three phase boundary lines meet, for each q case along 
with the ferromagnetic-paramagnetic phase boundary of the 
XY limit. 



tions is obtained by the convolution [24 
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where R({V(9ij)}) represents the bond moving and dec- 
imation given in Eqs.(3) and (4). For numerical practi- 
cality, the bond moving and decimation of Eqs.(3) and 
(4) are achieved by a sequence of pairwise combination 
of interactions, as shown in Fig. 1(c), each pairwise com- 
bination leading to an intermediate probability distribu- 
tion resulting from a pairwise convolution as in Eq.(5). 
We effect this procedure numerically, by generating 5,000 
interactions that embody the quenched probability dis- 
tribution resulting from each pairwise combination. The 
different thermodynamic phases of the model are iden- 
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FIG. 4: The calculated transition temperatures between the 
spin-glass phase and the paramagnetic phase, at p = 0.5, as 
a function of q, up to very large values of q = 720. Note the 
logarithmic scale of the horizontal axes. The slow decay of 
the transition temperature suggests that a zero-temperature 
spin-glass phase exists in the (q — s- oo) XY model limit. 



tified by the different asymptotic renormalization-group 
flows of the quenched probability distributions. Thus, 
we are able to calculate phase diagrams for any num- 
ber of clock states q. Our results are obtained by av- 
eraging over 30 to 50 different realizations of the initial 
±Jcos(9i — 9j) distribution into the 5,000 initial interac- 
tions. In this study, we consider even values of q and the 
calculated phase diagrams are symmetric around p = 0.5 
with the antiferromagnetic phase replacing the ferromag- 
netic phase, so that only the p — to 0.5 halves are shown 
below. 



III. CALCULATED PHASE DIAGRAMS FOR 
d = 3 q-STATE CLOCK AND XY SPIN GLASSES 

Our calculated phase diagrams for the q — 2,4,6,12 
clock spin- glass models are shown together in Fig. 2. 
The phase diagram for the XY limit, namely q — > oo, 
is also shown in Fig. 2, calculated here with q = 360 
clock states. In this limit, the spin-glass phase disap- 
pears at zero temperature, whereas the phase boundary 
between the ferromagnetic and paramagnetic phases nu- 
merically stabilizes, on the scale of the figure, for q > 6. 
The paramagnetic-ferromagnetic-spinglass reentrance as 
temperature is lowered, previously seen [25|,[26| for q = 2, 
namely the Ising case, is also seen here for the other q. 
As q is increased, it is found that the spin-glass phase 
retreats to lower temperatures while further protruding 
into the ferromagnetic phase. 

The calculated phase diagrams for the high-g mod- 
els, q = 12,18,36,360, are shown in Fig 3. As 
q is increased, the trend mentioned above, of the 
spin-glass phase retreating to lower temperatures while 
further protruding into the ferromagnetic phase, is 
also seen here. Furthermore, two new phenom- 
ena are simultaneously detected here: (1) dou- 



4 



0.05 



0.04 



0.03 



0.02 



■3 0.01 



0.04 



0.03 



0.02 



0.01 




-2 2 

Interaction U / <IUI> 



FIG. 5: (Color online) Asymptotic fixed distribution, under 
renormalization-group transformations, of the interactions in 
the spin-glass phase, namely the renormalization-group sink 
of the spin-glass phase. Note that these are strong-coupling 
distributions, as the average interaction strength < \U\ > 
diverges to infinity under the renormalization-group transfor- 
mations. Top: For the q = 6-state clock model in d — 3, 
for a trajectory starting in the spin-glass phase at p = 0.5 
and temperature 1/J = 0.05, the distributions after 20 and 
21 renormalization-group steps are shown. It is seen that 
these two distributions coincide, signifying a fixed distribu- 
tion. Bottom: Asymptotic fixed distributions of the spin- 
glass phases for the q = 2 (Ising), 6, 12-state clock mod- 
els in d = 3. These distributions are reached after 20 
renormalization-group steps, starting at p = 0.5 and tem- 
perature 1/J = 0.05. Note that the spin-glass sink fixed dis- 
tributions for different values of q coincide. The Lyapunov 
exponent is A = 1.93 for the single, universal distribution 
that is illustrated in this figure. 



ble reentrance, namely paramagnetic-ferromagnctic- 
spinglass-ferromagnetic phases as temperature is low- 
ered; (2) lateral reentrance, namely ferromagnetic- 
spinglass-ferromagnetic-paramagnetic phases as p is in- 
creased. Multiple reentrances have previously been seen 
in liquid crystal systems. [27l - l29| 

The (slow) disappearance of the g-state clock spin-glass 
phase is shown in Fig. 4, where the calculated spin-glass 
transition temperatures at p = 0.5 are shown as a func- 
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FIG. 6: (Color online) Comparison, showing the coinci- 
dence, of the chaotic visits of the consecutively renormal- 
ized interactions at a given position of the system (for 2000 
renormalization-group iterations) and of the asymptotic dis- 
tribution of the interactions across the system at a given 
renormalization-group step, for the spin-glass phase. The 
Lyapunov exponent is A = 1.93. Identical asymptotic be- 
havior occurs for all g-state clock spin-glass phases, as shown 
in Fig. 5. 



tion of q, up to very large values of q — 720. The slow 
decay of the transition temperature suggests that a zero- 
temperature spin- glass phase (30j exists in the q — > oo, 
namely XY model limit. 



IV. STABLE FIXED DISTRIBUTION AND 
CHAOTIC RENORMALIZATION-GROUP 
TRAJECTORY OF CLOCK SPIN-GLASS PHASES 

A. Stable Fixed Distribution 

For the spin-glass phase of the Ising model (q = 2), 
under repeated renormalization-group transformations, 
the quenched probability distribution of the interactions 
across the system becomes symmetric in ferromagnetic 
( Jij > 0) and antiferromagnetic ( Jij < 0) couplings, with 
the average magnitude of either type of interaction equal 
and diverging to infinity. [25] 

For the spin-glass phases of all g-state clock models, 
we find that under repeated renormalization-group trans- 
formations, the interaction values Vij(0) divide into two 
groups: 9 = 2irn/q and 8 = 2i:m/q, where n is an 
even integer, n=0,2,4,...,q-2, and m is an odd integer, 
m=l,3,5,...,q-l. The asymptotic renormalized quenched 
distribution of the interactions is symmetric, with inter- 
actions equal within each group mentioned above, but at 
each location overwhelmingly favoring one or the other of 
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FIG. 7: (Color online) The fixed distribution and, equiva- 
lently, chaotic renormalization-group trajectory onto which 
the phase boundary between the spin-glass and paramagnetic 
phases renormalizes, for the q = 2 and q = 6-state clock mod- 
els in d — 3. The Lyapunov exponent is A = 1.35. Identical 
asymptotic behavior occurs for all q. Note that the chaotic 
behavior and fixed distribution are at intermediate coupling 
strength. 



the two groups. Thus, the interaction difference between 
the two groups, 

j 9-2 1 9-1 

— V *A^n/q)- — %0W<7) = Uij, 

q ' Z n=0,2,... q ' Z ,n=l,3,... 

after many renormalization-group transformations, is 
randomly and equally distributed as positive or negative 
in our sampling of 5,000 interactions, which represent the 
distribution of interactions spatially across the system, 
with the average magnitude of either type of interaction 
equal and diverging to infinity. This asymptotic fixed 
distribution, namely the sink of the spin-glass phase, is 
shown in Fig. 5. Note that this asymptotic behavior is 
also consistent with the behavior of the Ising spin-glass 
phase (U reduces to J for the Ising case), recalled at 
the beginning of this section. This behavior leaves the 
system asymptotically frustrated. 



B. Chaotic Renormalization-Group Trajectory 

The fact that the Ising spin-glass phase is character- 
ized by the chaotic rescaling behavior of the interactions 
Q , and therefore of the correlations , is now well es- 
tablished P4l5| and is also seen here for the spin-glass 
phases of the g-state clock models in d = 3, as shown in 
Fig. 6. As with the Ising model Q , we have here a strong- 



FIG. 8: (Color online) Two different, non-coinciding strong- 
coupling fixed distributions: The fixed distribution and, 
equivalently, chaotic renormalization-group trajectory onto 
which the spinglass-ferromagnetic phase boundary of the 
q — 2-state model in d = 3 renormalizes. For comparison 
and distinction, the asymptotic fixed distribution and chaotic 
renormalization-group trajectory of the corresponding spin- 
glass phase is also shown. The latter curve is symmetric 
around Uij/ < \U\ >= 0, whereas the former curve is no- 
ticeably displaced towards positive (ferromagnetic) interac- 
tion values. Note that these are strong-coupling distributions 
and are therefore shown as a fraction of the diverging < \U\ >. 
Thus, even small shifts signify large differences in the interac- 
tion values. The Lyapunov exponents are A = 1.69 and 1.93, 
respectively for the phase boundary and phase sink cases. For 
q = 2, the form of the interaction in Eq.(l) is conserved under 
renormalization, which is reflected in the horizontal axis label 
here. 



coupling chaotic behavior: The values of the interaction 
difference Uij obtained by successive renormalization- 
group transformations at any specific location, divided 
by the average magnitude < \U\ > across the system, 
fall into a chaotic band. Thus, the Uij/ < \U\ > values 
are sampled within the band as shown in Fig. 6. The 
average magnitude < \U\ > diverges to infinity under 
repeated renormalization-group transformations. 

We thus realize that the spin-glass phases can be 
characterized by the Lyapunov exponent of the chaotic 
renormalization-group trajectory at any specific location, 



A = lim — In I — ^~ 
^ I dx 



k=0 



where x^ = Uij/ < \U\ > at step k of the renormaliza- 
tion-group trajectory. We have calculated the Lyapunov 
exponent A = 1.93 for the clock spin-glass phases of q = 
2, 6, 12 and presumably for the clock spin-glass phases of 
all q. 
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V. UNSTABLE FIXED DISTRIBUTIONS AND 
CHAOTIC RENORMALIZATION-GROUP 
TRAJECTORIES OF THE CLOCK 
SPINGL ASS-PARAMAGNETIC AND 
SPINGLASS-FERROMAGNETIC BOUNDARIES 

We find that the points on the various spin-glass phase 
boundaries also renormalize to a fixed distribution of 
the quenched interactions across the system and, equiv- 
alently, to a chaotic renormalization-group trajectory of 
the interaction at any single location in the lattice. The 
difference between the asymptotic rescaling behaviors in- 
side the spin-glass phase and on the spin-glass phase 
boundaries is that, under rescaling transformations, the 
fixed distribution and the chaotic trajectory are reached 
in a stable manner, with respect to initial conditions, for 
the spin-glass phase and are conversely unstable for the 
spin-glass phase boundaries. 
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FIG. 9: The fixed distribution and, equivalently, chaotic 
renormalization-group trajectory onto which the spinglass- 
ferromagnetic phase boundary of the q = 6-state clock model 
in d — 3 renormalizes. The coupled distributions of the inter- 
actions V(0)/ < \V(9)\ > are shown, while next to each axis, 
the fraction of points on the axis are given. The full fixed 
distribution has the form P({V(0), V(60), V(120), V(180)}), 
as a coupled function of its arguments. This is a strong- 
coupling behavior: The inset in the upper left panel shows, 
with the logarithmic vertical scale, the diverging < |V(#)| > 
as a function of renormalization-group iteration n, the con- 
secutive curves being for — 180, 120, 60, 0. It is seen that 
< \V{6)\ > diverges as b 0A6n for all 6. 



C. Equivalence of the Chaotic 
Renormalization-Group Trajectory and the 
Quenched Probability Fixed Distribution 



The distributions of the interaction difference Uij val- 
ues shown in Figs. 5 and 6, respectively obtained as 
the spatial distribution across the system after many 
renormalization-group transformations and the values 
obtained by successive renormalization-group transfor- 
mations at a specific location in the system, are in fact 
identical, as seen in Fig. 6. Thus, it is understood that 
the asymptotic fixed distribution is realized, after a given 
number of renormalization-group transformations, by the 
interactions at different locations being at different points 
of the same chaotic trajectory. 



A. The Spinglass-Paramagnetic Phase Boundary 

The phase boundary between the spin-glass and para- 
magnetic phases renormalizes to the fixed distribution 
and chaotic renormalization-group trajectory shown in 
Fig. 7. The interaction grouping, under rescaling, de- 
scribed before Eq.(6) also happens. However, this behav- 
ior here occurs at finite coupling < \U\ > = 0.686 for all q, 
in contrast to the asymptotic behaviors of the spin-glass 
phase (given above) and of the spinglass-ferromagnetic 
phase boundary (given below), which occur at strong 
coupling < \U\ >— > oo. 



B. The Spinglass-Ferromagnetic Phase Boundary 

The phase boundary between the spin-glass and fer- 
romagnetic phases renormalizes to a fixed distribution 
and chaotic renormalization-group trajectory at strong 
coupling < \U\ >— > oo. The interaction grouping, un- 
der rescaling, described before Eq.(6) does not happen. 
Thus, the interaction V(0) as a function of 6 has q/2 + 1 
different values. The asymptotic fixed distribution for 
q = 2 is shown in Fig. 8 and is characterized by the 
Lyapunov exponent A = 1.69. 

The asymptotic fixed distribution for q = 6 is shown 
in Fig. 9. The full fixed distribution has the form 
P({V(0), V(60), F(120), 7(180)}), as a coupled function 
of its arguments. The dominant configurations of this 
fixed distribution are shown in Table I. The system re- 
mains frustrated at all length scales. 



VI. CONCLUSION 

We have calculated, from renormalization-group the- 
ory, the phase diagrams of arbitrary even q-state clock 
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TABLE I: Dominant potentials in the asymptotic fixed 
distribution of the phase boundary between the spin-glass 
and ferromagnetic phases of the q — 6 clock model in d = 3. 
Thus, the system remains frustrated at all length scales. 



spin-glass models in d = 3. These models, for all non- 
infinite q, have a finite-temperature spin-glass phase, ex- 
hibiting a universal ordering behavior, independent of q. 
In addition to the chaotic rescaling behavior of the spin- 
glass phase, each of the two types of spin-glass phase 
boundaries displays, under renormalization-group trajec- 
tories, their own distinctive chaotic behavior, subdividing 
into two categories: strong-coupling chaos, in the spin- 
glass phase and distinctly on the spinglass-ferromagnetic 



phase boundary, and intermediate-coupling chaos, on the 
spinglass-paramagnetic phase boundary. We uniquely 
characterize each different phase and phase bound- 
ary exhibiting chaos by its distinct Lyapunov expo- 
nent, which we calculate. We demonstrate that, 
under renormalization-group, chaotic trajectories and 
fixed distributions are mechanistically and quantitatively 
equivalent. The spin- glass phases and the spinglass- 
paramagnetic phase boundaries exhibit universal fixed 
distributions, chaotic trajectories and Lyapunov expo- 
nents. In the XY model limit, our calculations indicate 
a zero-temperature spin-glass phase. 
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